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Abstract 

A series of tests were made in the NASA Lewis Icing 
Research Tunnel to determine how icing scaling results 
were affected by the choice of scale velocity. 
Reference tests were performed with a 53.3-cm-chord 
NACA 0012 airfoil model, while scale tests used a 
27.7-cm-chord 0012 model. Tests were made with 
rime, mixed and glaze ice. Reference test conditions 
included airspeeds of 67 and 89 m/s, an MFD of 40 pm, 
and LWCs of 0.5 and 0.6 g/m 3 . Scale test conditions 
were established by the modified Ruff (AEDC) scaling 
method with the scale velocity determined in 5 ways. 
The resulting scale velocities ranged from 85 to 220% 
of the reference velocity. This paper presents the ice 
shapes that resulted from those scale tests and compares 
them to the reference shapes. It was concluded that for 
freezing fractions greater than 0.8 as well as for a 
freezing fraction of 0.3, the value of the scale velocity 
had no effect on how well the scale ice shape simulated 
the reference shape. For freezing fractions of 0.5 and 
0.7, the simulation of the reference shape appeared to 
improve as the scale velocity increased. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Model chord, m 

hfiim Water-film thickness at leading edge, m 

K 0 Modified inertia parameter, dimensionless 

LWC Cloud liquid-water content, g/m 3 

M Mach number, dimensionless 

MVD Water droplet median volume diameter, pm 

n Freezing fraction, dimensionless 

p lo , Total pressure, nt/m 2 

*Senior Research Associate, Member AIAA 
1 Associate Professor, Associate Fellow AIAA 


p s , Static pressure, nt/m 2 

Re Reynolds number of model, dimensionless 

Re g Reynolds number of water droplet, 

dimensionless 

t s , Static temperature, °C 

t, ot Total temperature, °C 

V Air velocity, m/s 

We Weber number based on droplet size and water 
properties, dimensionless 

4> Droplet energy transfer tenns in energy 

balance, °C 

p. Ice density, g/m ’ 

9 Air energy transfer terms in energy balance, 

C 

T Accretion time, min 

General Subscripts 

r Reference condition 

5 Scale condition 

Introduction 

This paper presents the results of a series of tests 
recently made in the NASA Lewis Icing Research 
Tunnel to determine how scaled ice shapes were 
affected by the choice of scale velocity. Scale velocity 
was found by applying four physically based methods. 
A fifth velocity was simply an average of two of the 
others. One physical method matched the scale and 
reference values of water-film thickness, one matched 
the velocity, one, the Weber number and one, the 
Reynolds number. By establishing which of these 
approaches was most effective, it was hoped to gain 
insight into the importance of the physics on which the 
method was based. 

A number of scaling methods have been developed to 
address the problem of testing models in facilities with 
limitations in either model size or test conditions. 
These methods have all been based on a physical 
description of the ice accretion process with various 
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assumptions and simplifications. The usual approach in 
developing a scaling method is to write expressions for 
the water droplet collection efficiency, total ice 
accretion on the surface and an energy balance at the 
stagnation line of the airfoil. The scale and reference 
values of these expressions are equated to solve for the 
scale test conditions. 

Langmuir and Blodgett 1 defined the modified inertia 
parameter, K 0 , which is directly related to the local 
collection efficiency at the stagnation line. This 
dimensionless parameter has been widely used in 
scaling methods to insure the scale droplet trajectories, 
and, thus, collection efficiencies, match the reference. 
The non-dimensional ice accretion is represented by the 
accumulation parameter, A c = (LWC V x)/(c p,). Heat 
transfer characteristics can be included by writing the 
Messinger energy balance 2 in terms of the freezing 
fraction, n, the water energy transfer parameter, p, and 
the air energy transfer parameter, 9. (f> and 9 are 
collections of terms in the energy equation that relate to 
the transfer of energy due to impinging droplets and to 
airflow over the model, respectively. Ruff 1 perfonned 
scale tests with several combinations of these 
parameters and achieved the best scaling results when 
all five were matched to the reference values. This 
approach to scaling is known as the Ruff (or AEDC) 
Method. 

By matching five scaling parameters, the Ruff Method 
permits the calculation of five scale test conditions from 
the corresponding reference conditions. They are: 
static temperature, t sl , cloud water droplet median 
volume diameter, MVD, cloud liquid-water content, 
LWC, spray duration, X, and test-section total pressure, 
p m . For tunnels, such as the IRT, which do not control 
the test-section pressure, a modified Ruff Method has 
been used 4 in which the air-energy transfer parameter, 
9, is ignored, and the test-section pressure is determined 
by ambient conditions and test-section velocity. 

In the Ruff Method, as with several other scaling 
methods, the user is free to specify the scale velocity. 
One way to specify the scale velocity is to equate the 
scale to the reference velocity. Besides being 
convenient, this approach insures that 9 and the Mach 
number, M, match the reference values. Other ways to 
choose the scale velocity can be devised by considering 
additional physical processes and parameters that may 
be important in determining ice shape. 

In 1986, Olsen and Walker 5 performed tests with close- 
up photography of the icing surface. This study 
provided evidence that physics of the surface water 
might be important in addition to heat transfer in 
determining the shape of glaze ice. About the same 


time, Bilanin 6 published a study of scaling in which he 
argued that a number of parameters, including the 
Weber number, We, and Reynolds number, Re, need to 
be considered in formulating scaling methods. Because 
We is the ratio of aerodynamic to surface tension forces, 
it was believed to represent at least some of the surface- 
water effects and was chosen by Bilanin and 
Anderson 7 ’ 8 to be included in a new scaling method. 
This constant- We method is in fact the Ruff method 
with the added restraint that the scale and reference We 
match, and from this requirement, the scale velocity can 
be found. For subscale tests, the constant- We method 
produces a scale velocity higher than the reference. 

A third method is to match the scale and reference Re, 
following Bilanin. 6 Because the Reynolds number 
controls a number of physical processes, including both 
flow and heat transfer, it is logical that the scale value 
should be set equal to the reference value. The use of 
the constant-/?e approach in subscale tests produces a 
scale velocity still higher than that found by matching 
We. Although very high scale velocities are rarely 
practical, the method needs to be evaluated to 
determine the importance of Re to ice-accretion scaling. 

A method that gives a lower-than-reference scale 
velocity for subscale models was proposed by Feo. 9 He 
suggested matching the scale and reference leading- 
edge water-film thickness, normalized by model size. 
Feo’s water-film thickness was based on a correlation 
of data taken as part of a study of the effects of rain on 

°c LWC 1/2 (Re)~ ll4 (Re 8 ) 15/4 (l Ve)~ 9 ' 4 
c 

aircraft performance. 10 The best fit was found to be 

This correlation was developed from tests with larger 
water droplets and higher liquid-water contents than 
those experienced in normal icing encounters. 
However, in the absence of more relevant data, it can be 
taken as an approximate representation of the thickness 
of the water film which forms during glaze ice 
accretion. The ratio of the scale airspeed to the 
reference can be found by equating the scale and 
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reference non-dimensional water film thickness: 

Here it was assumed that the droplet velocity which 
appears in We and Reg is proportional to the airspeed, 
and that differences between scale and reference static 
pressure and water properties are negligible. 
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Figure 1 shows typical scale velocities found by the 
four methods described above for a reference velocity 
of 67 m/s and for a !4-size scale model. A fifth scale 
velocity is also shown; it is the average of the constant- 



ly We V Re 

Method 

Figure 1. Scale Velocities for W scale Model. 

Reference Velocity, 67 m/s. 

We and constant-// c> velocities. This average- V method 
has no physical basis, but it can be used as part of a 
sequence of tests to assist in spotting possible trends 
with velocity. 

In the present study, scaled tests with half-size models 
were made using the Ruff Method with scale velocity 
found by each of these five methods. 

Reference Conditions 

Five reference cases with full-size models were chosen 
for study. These are shown in table I. Case 73 


Table I. Reference Test Conditions 
Reference Model: 53.3-cm-Chord NACA 0012 


Case 

t St) 

°c 

V, 

m/s 

MVD, 

pm 

LWC, 

g/m 3 

T , 
min 

n 

73 

-20.6 

89 

40 

0.5 

10.0 

1.0 

110 

-13.9 

67 

40 

0.6 

11.2 

0.8 

111 

-11.1 

67 

40 

0.6 

11.2 

0.7 

112 

-8.4 

67 

40 

0.6 

11.2 

0.5 

113 

-5.6 

67 

40 

0.6 

11.2 

0.3 


produced rime ice, while 1 10-1 13 gave mixed or glaze. 
Cases 110-113 had the same velocity, spray conditions 
and time, with the temperature increasing from case 
1 10 to case 1 13 to provide a range of freezing fractions. 
The reference MVD and LWC were selected so that 
both reference and scale conditions could be tested with 
the same nozzle set in the IRT. Spray times were 
chosen to provide the same accumulation parameter for 
all cases. 


Test Description 

NASA Lewis Icing Research Tunnel 

Figure 2 is a plan view of the NASA Lewis Icing 
Research Tunnel (IRT). 11 The IRT has a test section 



Figure 2. NASA Lewis Icing Research Tunnel (IRT). 

width of 2.7 m and a height of 1.8 m. It has a 
refrigeration system that allows accurate control of the 
test-section temperature from -40 to 4°C. A water 
spray system with ten spray bars simulates the 
conditions in a natural icing cloud. The test-section 
cloud droplet size, MVD, and liquid-water content, 
LWC, depend on spray-bar air and water pressures. The 
relationships among these pressures, the tunnel airspeed 
and the cloud properties are established periodically by 
a series of tunnel calibration tests. 12 The cloud has 
been calibrated over a range of test-section airspeeds 
from 22 to 156 m/s and droplet median volume 
diameters from 14 to 50 pm. Two sets of spray 
nozzles, the Mod-1 and Standard nozzles, are used to 
provide different ranges of liquid-water content. 
Depending on the nozzle set, the airspeed and the 
droplet diameter, the test-section liquid-water content 
can be controlled from less than 0.2 to over 5 g/m 3 . 
The Mod-1 nozzles were used for all tests reported 
here. 

Test Models and Procedures 

Ice accretion was measured on NACA 0012 airfoils 
with chords of 53.3 cm (reference model) and 26.7 cm 
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(scale model). Each was mounted vertically, at 0° 
angle of attack, in the center of the IRT test section. 
The models were machined from solid aluminum. Both 
had a span of 61 cm and were placed between end 
plates as shown in figure 3. The models’ mid spans 



Figure 3. 53.3-cm-Chord NACA 0012 Model in 
IRT Test Section. 

were located on the tunnel centerline (midway between 
floor and ceiling of the test section). Five horizontal 
lines were marked around the leading edge of each 
airfoil to indicate the location of the mid-span and of 
positions 2.5 and 5 cm above mid-span and 2.5 and 5 
cm below mid-span. Ice shapes were recorded at all but 
the 5-cm-below-mid-span location. All shapes reported 
were recorded by hand tracing. 

The IRT spray system recirculates water to permit the 
stabilization of air and water pressures prior to opening 
the water valve at each nozzle. This feature virtually 
eliminates start-up transients that were present 
previously in the IRT. Once the desired temperature 
and airspeed had been achieved and spray-bar 
conditions had been set and stabilized, the water valves 
were opened to initiate the spray, and the spray timer 
was started. At the completion of the icing spray time, 
the spray was turned off and the tunnel fan was 
stopped. Personnel then entered the test section, and a 
thin heated plate was used to cut horizontal slices into 
the accreted ice at the desired span-wise locations. The 
shapes were traced onto cardboard templates and later 


digitized for computer storage. Finally, the model was 
cleaned and the procedure repeated for the next test. 

Average Test Conditions 

Actual reference and scale conditions for the tests 
reported here are given in table II along with the 
corresponding scaling parameters. Tunnel and cloud 
conditions were recorded at 1-sec intervals for tests 
with the 26.7-cm-chord model and at 2-sec intervals for 
the 53.3-cm-chord model. Each test condition was also 
determined from more than one instrument. The total 
temperature was recorded with eleven type T 
thermocouples located on the turning vanes just 
downstream of the heat exchanger. Two pitot-static 
probes on opposite walls at the entrance to the IRT test 
section were used to determine test-section velocity. 
The MVD and L WC were calculated from the spray-bar 
air pressure and water-air differential pressure for each 
of the ten spray bars. The test conditions in Table II 
have been averaged over the multiple instrument 
readings as well as over the duration of the spray. 

Uncertainty in Parameter Values 

Typically, the total temperature drifted up and down 
during each test by about ±0.1 °C. Fluctuations in the 
velocity and the LWC on the order of ±1% were 
observed, while the droplet size showed little variation 
after the first 20 sec or so. 

Estimates of the uncertainty in the reported conditions 
were made by considering fluctuations of the values 
with time, possible instrument errors including 
calibration, uncertainties in tunnel calibration of MVD 
and LWC, and differences in measurements from one 
location to another. It was concluded that for the 
conditions of the tests reported here total temperatures 
are probably known to within ±1.5 °C, velocity, ±2.5%, 
MVD, ±11% and LWC, ±1 2%. 

Results and Discussion 

The ice shapes recorded in these tests will be shown 
with coordinate axes graduated in cm. These axes 
represent the full-size coordinates of the ice shapes for 
the reference model. The coordinates for the '/ 2 -scale 
shapes were doubled to pennit direct comparison with 
the reference. 

Repeatability tests were made to establish typical run- 
to-run differences in ice shape. Checks of the span- 
wise ice-shape uniformity were also made, because a 
highly non-uniform cloud will produce ice shapes 
which are very sensitive to tracing location. 

Repeatability of Ice Shapes 

Figure 4 presents three comparisons of ice formed on 
different occasions at the same test conditions and with 


NASA/CR— 2003-21 1827 


4 



-2 0246 -2 0246 -2 0246 

(a) 3-26-98 Run 2 (b) 3-26-98 Run 2 (c) 3-26-98 Run 2 

- 10-26-95 Run 3 12-11-95 Run 2 3-27-98 Run 2 

Figure 4. Repeatability of Ice Shapes, 1995-1998. Case 1 13 Reference at Mid-Span. 53.3-cm-Chord NACA 0012 
Airfoil. t st , -3°C; V, 67 m/s; MVD, 40 pm; LWC, .65 g/m 3 ; x, 1 1.2 min. Axis Units, cm. 


the same 53.3-cm-chord NACA 0012 model. All ice 
shapes were recorded at the vertical center of the IRT 
test section. Figure 4(a) compares shapes recorded on 
3-26-98 and 10-26-95. The two agree very closely 
except for some minor differences in the lower horn. 
This agreement is impressive considering that the 
tunnel spray bar system had been replaced and 
recalibrated in 1997. A second comparison is given in 
figure 4(b) in which the 3-26-98 shape is compared 
with one from 12-11-95. The latter shape is 
significantly smaller, and the impingement limits are 
nearer the leading edge. The final comparison, in 
figure 4(c), is between the 3-26-98 test and one made 
the following evening. Again, excellent repeatability of 
the ice shape is evident. It is clear from these 
comparisons that the IRT is capable of excellent 
repeatability even between tests performed years apart. 
Occasionally, however, test results can vary from one 


test entry to another. All scaling results reported here 
were made during the same tunnel entry to minimize 
variations of ambient conditions and tunnel 
characteristics that might influence ice shape. 


Span-wise Unifonnit 


Ice shapes recorded at different positions along the 
model span are shown in figure 5. The mid-span 
position is at the vertical center of the tunnel test 
section. Figure 5(a) compares the mid-span shape with 
the tracing taken 5 cm above the mid-span position. 
Figure 5(b) compares the mid-span and 2.5 cm above 
mid-span, and figure 5(c) shows the mid-span shape 
and one taken 2.5 cm below mid-span. The excellent 
agreement among the tracings taken at the four 
positions demonstrates both the consistency of the 
tracing process and the uniformity of the IRT cloud 
over at least this limited vertical range. The uniformity 



(a) Mid Span (b) 

5 cm Above Mid Span 


Mid Span (c) 

2.5 cm Above Mid Span 


Mid Span 

2.5 cm Below Mid Span 


Figure 5. Span-wise Uniformity of Ice Shape. Case 1 13 Reference, 3-26-98 Run 2. 53.3-cm-Chord NACA 0012 
Airfoil, t st , -3 °C; V, 67 m/s; MVD, 40 pm; LWC, .65 g/m 3 ; x, 1 1.2 min. Axis Units, cm. 
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indicated by figure 5 was typical of 
both scale and reference tests for 
this study. These results gave 
confidence that scale and reference 
shape comparisons would not be 
affected by small errors in tracing 
location. Scale-to-reference 

differences in shape reported in the 
past may sometimes have been due 
in part to non uniformity in the 
cloud such that tracings taken at 
slightly different span-wise 
locations might have shown better 
agreement. 



(a) Feo; V s , 76 m/s 



(b) Constant- Velocity; V s , 89 m/s. 


(c) Constant-Weber-Number; V„ 
118 m/s. 


Rime Ice (Case 73) 

For rime ice water freezes 
immediately on impact because 
convective heat transfer absorbs all 
the latent heat released on freezing. 

Consequently, a description of the 
heat balance at the surface is 
unnecessary, and, with no liquid 
water on the ice, dynamics of a 
liquid surface are irrelevant. Thus, 
only the accumulation and inertia 
parameters, A c and K 0 , have to be 
matched between scale and 
reference, and all of the methods to 
choose scale velocity should have 
been equally successful in 
simulating the reference ice shape. 

An inability of scaled shapes to 
simulate the reference shape is 
therefore evidence of a problem with either the MVD 
(affecting Kg ) or LWC (affecting A c ) calibration of the 
cloud. 

Figure 6 gives the reference and scale ice shapes for the 
rime case 73 tests for the Feo (fig. 6(a)), constant- V 
(fig. 6(b)), constant- We (fig. 6(c)) and average- V 
(fig. 6(d)) methods of selecting scale velocity. The 
scale velocity required to satisfy constant-//^ was 
219 m/s, well in excess of the IRT capability, so this 
method could not be tested for this case. Both the Feo 
and constant- V methods produced scale shapes that 
matched the reference in shape and quantity. Flowever, 
while the constant- We and average- V results simulated 
the reference shape, neither produced sufficient 
quantity of ice. It was suspected, therefore, that the 
MVDs for the results shown in figures 6(c) and (d) were 
properly scaled, but the scale LWC values were too low. 




(d) Average-Velocity; V s , 156 m/s. 


Figure 6. Scaling Results for Rime Ice Using Four Methods of Choosing 
Scale Velocity. Case 73 (n = 1.0). Test Conditions in table II. Axes 
Units in cm. Scale Ice Shape Coordinates Doubled. 

Reference Shape, c r , 53.3 cm; V r , 90 m/s. 

Scale Shape; c s , 27.7 cm; V s , as indicated. 

of an annual check. Slightly different procedures were 
used than in past calibrations, and the LWCs for given 
spray-bar pressures were found to be somewhat lower 
than for the 1997 calibrations. The greatest difference 
in LWC between the two calibrations was found at the 
highest velocities. When this new calibration was 
applied to the March data, the leading-edge thicknesses 
of these rime shapes correlated with the new 
accumulation parameter as shown in figure 7. Figure 7 
includes data from the rime tests reported in figure 6 as 
well as from additional rime shapes that will not be 
presented here. Because the 1998 calibration appeared 
to explain the observed results with rime ice, it was 
used to determine the LWCs for the March 1998 tests as 
given in table II. This experience showed that a series 
of scaling tests perfonned with rime ice conditions can 
be an effective tool to check the consistency of icing 
tunnel cloud calibrations over a range of velocity. 


The target MVD and LWC values for these March 1998 
tests were established by setting spray-bar air and water 
pressures according to a 1997 calibration of the IRT. 
The tunnel was recalibrated in November 1998 12 as part 


The Mach number for the average- V scaled conditions 
(fig. 6(d)) approached 0.5, well in excess of the 
compressible limit 13 of about 0.3. However, figure 6(d) 
does not indicate a scaled ice shape any different from 
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Figure 7. Correlation of Maximum Rime Ice Thickness 
with Accumulation Parameter. A c Based on 
1998 IRT Cloud Calibration. 

that in figure 6(c) for which M= 0.37. This result and 
those to be given later suggest that compressibility 
probably has little effect on ice accretion, at least for 
M< .5. 

Mixed Ice (Cases 1 10 and 111) 

The case 110 (n = 0.8) reference ice shape was mostly 
rime but with a narrow glaze region at the leading edge, 
as can be seen in figure 8. Due to time constraints, 
scaling tests for case 110 were made only with scale 
velocities from the Feo, constant- V and constant- We 
methods. Results are shown in figures 8(a), (b) and (c). 
For this icing condition, the quality of the scaled shape 
relative to the reference appeared to be independent of 
which of these three methods were used. All three 
scaled shapes, when coordinates were doubled, were 


approximately the same size as the reference, but none 
exactly reproduced the leading-edge features. 

Figure 9 presents the scaling results for case 111 (/? 
= 0.7). The reference shape for this case featured ice 
with small horns. A large feather structure immediately 
adjacent to the main shape grew into the side of the 
horns. For this case, scaling tests were only made with 
the constant- We, average- V and constant -Re methods. 
None of the scaled ice shapes matched all the 
prominent features of the reference shape. The 
constant- We and average- V methods (figs. 9(a) and (b)) 
failed to reproduce the main horns although the large 
adjacent feathers were approximated. 

The constant-Re-scaled ice results are shown in 
figure 9(c). Although the scale Mach number for this 
test was 0.47, well above the maximum for 
incompressible flow, the scale ice shape simulated the 
reference shape in the region of the leading edge. 
However, the abundant feather growth of the reference 
ice was nearly non existent in the scaled ice. It is 
possible that the high scale velocity prevented 
significant feather growth by causing fragile structures 
to shed. 

Glaze Ice (Cases 112 and 113) 

Case 112 had a freezing fraction of 0.5. The reference 
ice shape featured the horns typical of glaze ice with 
relatively small feathers farther back on the airfoil. All 
five scale-velocity methods were tested. 

The Feo method (fig. 10(a)) produced scaled ice with a 
convex shape at the leading edge region. The resulting 
horns were somewhat farther back than the reference 
horns although the hom size was simulated well. As 
the scale velocity was increased by applying the 
constant-R (fig. 10(b)) and constant-IRe (fig. 10(c)) 
methods the horns moved farther forward. The lower 



(a) Feo; V s , 57 m/s. 



(b) Constant-Velocity; V s , 67 m/s. 



(c) Constant-Weber-Number; V s , 
88 m/s. 


Figure 8. Scaling Results for Mixed Ice Using Three Methods of Choosing Scale Velocity. Case 1 10 (n = 0.8). Test 
Conditions in table II. Axes Units in cm. Scale Ice Shape Coordinates Doubled. 

Reference Shape; c r , 53.3 cm; V r , 67 m/s. 

Scale Shape; c s , 27.7 cm; V s , as indicated. 
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(a) Constant-Weber-Number; V s , (b) Average-Velocity; V s , (c) Constant-Reynolds-Number; 

88 m/s. 117 m/s. V s , 1 47 m/s. 

Figure 9. Scaling Results for Mixed Ice Using Three Methods of Choosing Scale Velocity. Case 111 (n = .7). Test 
Conditions in table II. Axes Units in cm. Scale Ice Shape Coordinates Doubled. 

Reference Shape, c r , 53.3 cm; V r , 67 m/s. 

Scale Shape, c s , 27.7 cm; V s , as indicated. 



(a) Feo; V s , 57 m/s. 


(b) Constant-Velocity; V s , 67 m/s. (c) Constant-Weber-Number; V s , 

87 m/s. 



-2 0246 -2 0246 


(d) Average-Velocity; V s , 117 m/s. (e) Constant-Reynolds-Number; 

K, 147 m/s. 

Figure 10. Scaling Results for Glaze Ice Using Five Methods of Choosing 
Scale Velocity. Case 112 (/? = 0.5). Test Conditions in table II. 
Axes Units in cm. Scale Ice Shape Coordinates Doubled. 

Reference Shape, c r , 53.3 cm; V r , 67 m/s. 

Scale Shape, c s , 27.7 cm; V s , as indicated. 


velocity to that of the constant- V 
method (fig. 10(d)), the scale test 
produced a good simulation of all 
features of the reference shape. 

The final increase in scale 
velocity, to that of the constant- 
Re method, produced the scale 
ice shape of figure 10(e). This 
shape had a deep trough between 
the horns which did not exist in 
the reference shape, and the 
horns extended into the flow 
farther than for the reference. 
Table II shows that the scale 
Mach number was in the 
compressible regime, but results 
cited previously suggested that M 
should have little significant 
effect on the leading-edge ice 


horn of the constant- We shape simulated the reference 
well, but the upper horn location was still somewhat aft 
of the reference. With a further increase in scale 


shape. The scale freezing 
fraction was about 10% lower than the reference, but 
probably more important, the total temperature was just 
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above the freezing temperature of water. This high 
temperature is probably the cause of the leading-edge 
distortion of the scaled shape. 

In case 1 13 (fig. 1 1), the freezing fraction was 0.3. The 
reference ice shape was convex except for a small dip at 
the leading edge. The upper and lower horns grew 
nearly normal to the undisturbed airflow. Scaling 
results were obtained for the Feo, constant- V, constant- 
ly and average- V methods. The constant-Re method 


of case 1 12 (fig. 10(e)). Note from table II that the total 
temperatures of these two scale tests were nearly equal, 
a few tenths of a degree above freezing. The similarity 
of the features of these shapes along with the 
coincidentally high total temperatures support the 
supposition that these ice shapes were heavily 
influenced by the temperature. 

Concluding Remarks 




(a) Feo; V s , 57 m/s. 


(b) Constant Velocity; V s , 67 m/s. 



(c) Constant Weber Number; V s , 
88 m/s. 



(d) Average Velocity; V s , 118 m/s. 


Figure 11. Scaling Results for Glaze Ice Using Four Methods of Choosing 
Scale Velocity. Case 113 (n = .3). Test Conditions in table II. 
Axes Units in cm. Scale Ice Shape Coordinates Doubled. 

Reference Shape, c r , 53.3 cm; V r , 67 m/s. 

Scale Shape, c s , 27.7 cm; V s , as indicated. 


gave a scale total temperature of 3.2 °C; therefore, it 
was not tested. The scale shapes for the Feo (fig. 1 1(a)) 
and constant- V (fig. 11(b)) methods were similar; both 
approximated the leading-edge region of the reference 
ice shape. However, the upper horns of the scale 
shapes were located too far back and the lower horns 
were smaller than the reference. The scaled shape for 
the constant- We approach (fig. 11(c)) simulated the size 
and location of the upper horn well, but not the lower 
portion of the leading-edge region. The average- V 
method (fig. 11(d)) produced a scaled shape with 
characteristics much like those of the constant-//^ test 


To evaluate ways to choose scale 
velocity, a series of scaling tests 
were performed using the modified 
Ruff scaling method with 5 
approaches to selecting scale 
velocity. Reference tests used a 
53.3-cm-chord NACA 0012 
model, and scale tests, a 27.7-cm- 
chord NACA 0012. Scale 
velocities ranged in value from 
85% of the reference velocity 
(using a Feo constant- water-film 
thickness) to 220% (matching the 
scale Re to the reference.) The 
general conclusions from these 
tests can be summarized as 
follows: 

A series of rime scaling tests 
using a range of scale 
velocities can be used to 
identify inconsistencies in 
tunnel cloud calibration. 

For rime ice and ice formed 
with freezing fractions as low 
as 0.8, there appeared to be 
little effect of scale velocity on 
ice shape. 

For ice formed with a freezing 
fraction of 0.7, limited 
evidence suggested that 
improved simulation of the 
reference shape at the leading 
edge may have been achieved as the scale velocity 
increased. However, for the scale tests using the 
constant-/^ approach, the feather structure was 
inadequate to simulate that of the reference ice 
shape. The high velocity of this test may have 
hindered feather development. 

4. For ice formed with a freezing fraction of 0.5, 
improved simulation of the reference shape was 
achieved as the scale velocity increased, as long as 
scale total temperature did not exceed 0 °C. 
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5. At the lowest freezing fraction tested, n = 0.3, the 
scale velocity appeared to have little effect on the 
scaled ice shape, but few cases could be run 
without exceeding a total temperature of 0 °C. 

6. Tests in which the scale total temperature exceeded 
0 °C produced scaled ice shapes whose leading- 
edge features displayed no similarity to those of the 
reference shape, which had been formed at a lower 
temperature. 

7. Mach numbers as high as 0.5 appeared to have no 
effect on the main ice shape. 

The conclusions presented here are tentative. To insure 
the validity of the ice shapes on which they are based, 
the tests need to be repeated utilizing the most recent 
IRT cloud calibration. 

Neither the Feo constant-water-film-thickness nor the 
constant- We approach to choosing the scale velocity 
appeared to provide consistently better scaling than 
other methods tested. This result suggests that if the 
physics of surface water films is important to 
determining glaze ice shapes, neither of these methods 
adequately captures that physics. 

Some of the results of this study may point to the 
importance of Re in scaling. Re has generally been 
neglected because attempts to match the scale and 
reference Re often result in scale velocities or 
temperatures that are excessive. Practical scaling 
methods often require compromises, but these can best 
be applied only after the most important parameters 
have been identified. Further evaluation of constant-Re 
scaling is therefore encouraged. 

Further study of the heat-transfer terms n. (]), and 9 may 
also be warranted, with the objective of detennining if 
the expressions for these parameters properly describe 
the physics of ice accretion. 
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